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Abstract. The Pfaff lattice is an integrable system arising from the SR-group 
factorization in an analogous way to how the Toda lattice arises from the QR-group 
factorization. In our recent paper [Intern. Math. Res. Notices, (2007) rnml20], we 
studied the Pfaff lattice hierarchy for the case where the Lax matrix is defined to be 
a lower Hessenberg matrix. In this paper we deal with the case of a symplectic lower 
Hessenberg Lax matrix, this forces the Lax matrix to take a tridiagonal shape. We 
then show that the odd members of the Pfaff lattice hierarchy are trivial, while the 
even members are equivalent to the indefinite Toda lattice hierarchy defined in [Y. 
Kodama and J. Ye, Physica D, 91 (1996) 321-339]. This is analogous to the case of 
the Toda lattice hierarchy in the relation to the Kac-van Moerbeke system. In the 
case with initial matrix having only real or imaginary eigenvalues, the fixed points of 
the even flows are given by 2 x 2 block diagonal matrices with zero diagonals. We also 
consider a family of skew-orthogonal polynomials with symplectic recursion relation 
related to the Pfaff lattice, and find that they are succinctly expressed in terms of 
orthogonal polynomials appearing in the indefinite Toda lattice. 
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1. Introduction 



The Toda lattice hierarchy of nonlinear differential equations may be formulated as the 
continuous limit of the QR-algorithm for diagonalizing a symmetric matrix [T9J, [8j 20J. 
There exists a number of generalizations of the QR-algorithm for diagonalizing matrices 
while preserving different symmetries, the HR and SR-algorithms are two that we will 
consider in this paper [H [7J [TOj H9) . In particular the SR-algorithm relies on SR- 
factorizations whose goal is to factor an invertible matrix into a symplectic matrix and 
an upper triangular matrix [6j. The Pfaff lattice hierarchy was recently introduced to 
describe the partition functions of the orthogonal and symplectic ensembles of random 
matrices, and equivalently, to give the evolution of skew-orthogonal polynomials (see 
[U, [3]). The Pfaff lattice hierarchy may be viewed as a continuous limit of the SR- 
algorithm for diagonalizing a lower Hessenberg matrix, and in fact was originally 
discovered in this context [22] . We showed in our recent work [16] that the Pfaff lattice 
is an integrable system in the Arnold-Liouville sense, and in the case that the initial 
conditions are not lower Hessenberg matrices from the symplectic algebra, sp(n), we 
showed that the fixed points of the Pfaff flows with an initial condition having real 
distinct eigenvalues are given by 2 x 2 block diagonal matrices with zero diagonals. 

In this paper we will consider the Pfaff lattice flow with an initial condition given 
by a matrix from the symplectic algebra sp(n). The Lie algebra sp(n) is the set of 



2n x 2n matrices satisfying JL + L T J = with 

j = di ^((_° 1 i))= 7 »®(-i I 

Here diag2(A) denotes the 2x2 block diagonal matrix of A, i.e. diag 2 (A) = 
diag 2 (Ai, . . . , A 2 ) with the 2x2 blocks Aj in the diagonal of A, and I n is the n x n 
identity matrix. 

Remark 1.1 Matrices in the symplectic algebra, sp(n), are called Hamiltonian in the 
linear algebra and control theory literature, we will avoid this nomenclature to avoid 
confusion with the Hamiltonians which are constants of motion and appear throughout 
the paper. The term symplectic algebra is also the one commonly used in the literature 
on the Pfaff lattice and in representation theory J77] /. 

The spectrum of a real symplectic matrix comes in three types: 

(a) Pairs of real eigenvalues (z, —z). 

(b) Paris of imaginary eigenvalues (z, —z). 

(c) Or quadruples of complex eigenvalues ( Z • Z ^ Z j Z J • 

Symplectic matrices appear in a number of control theory problems most notably in 
solving the real algebraic Riccati equation (see e.g. [9]). Our J is a permutation of the 
traditional J, say J, used in the literature on the Hamiltonian eigenvalue problem, i.e. 
1 
-1 

is defined on lower Hessenberg matrices (see below). 

The Pfaff lattice can be viewed as an sp-version of the Toda lattice, it has the 
following form with sp-projection [H [16], 



J 



In ■ We use this J because of its connection to the Pfaff lattice which 



81 

— = [^(VH),L}. 

The Pfaff lattice is associated with the Lie-Poisson structure induced by the splitting 
sl(2n) = t(Bsp(n). Recall that the splitting gives the following projections of an element 
X E *[(2n,R), 



7r e (X) = X_ - J (X + f J+-(X - JX T J) 



7T 



(X) = X + + J (X + ) T J + ^ (X + JX'^J) 



Bp 



;i.2) 



where X is the 2x2 block diagonal part of X, X + (resp. X_) is the 2x2 block upper 
(resp. lower) triangular part of X. In particular, we have 

( ( hj 2 2 ... 2 \ 

-k h 2 h • • • 2 
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k=l 
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where I2 is the 2x2 identity matrix, and O2 is the 2x2 zero matrix. 

As in the case of the Toda lattice hierarchy, one can also define the Pfaff lattice 
hierarchy (see [16], and also Section [3]), 
dL 

[7T SP (V^),L] for j = l,...,2n-l, (1.4) 



where Hj 



1 



dti 



and L is a Hessenberg matrix given by, for example with n — 3, 
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The Cfe are invariants of the Pfaff lattice, in the traditional Pfaff lattice consider by 
Adler et. al. [1] the Ck are chosen to be 1. It was shown in [16] that the c& are 
Casimirs of the Pfaff lattice. In [16] it was shown that for L with distinct eigenvalues 
the Pfaff lattice hierarchy is an integrable system in the Arnold-Liouville sense. In 
addition it was shown that if L is not symplectic and has distinct real eigenvalues, then 
as tj — > 00 for any j, there is a diagonal matrix P(tj) such that the normalized matrix 
L(tj) := P(tj)L{tj)P(tj)~ l approaches a 2 x 2 block upper triangular matrix, where the 
diagonal blocks are sorted by the size of z'-> for the real eigenvalues z appearing in each. 
In this paper, we deal with the Pfaff lattice when the initial matrix is symplectic, and 
discuss the dynamical structure in a connection with the Toda lattice. 

A natural question at this point is then what happens when we restrict the Pfaff 
lattice variables to matrices which are both symplectic and lower Hessenberg. The first 
consequence is that such matrices have the compact form: 

/ ci \ 
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that is, L has a 2 x 2 block tridiagonal form L 
having hj = O2 for \i — j\ > 1 and 



1.6) 



j )l<i,j<n 



J 



with 2x2 block matrices 



k,k 




k,k+l 



I 



k+l,k 




where Ck = ±1. We call a matrix in this form iJ 5 -tridiagonal. This is analogous 
to restricting the Toda lattice variables to a matrix which is tridiagonal and skew- 
symmetric (i.e. in so(n)); a situation we will remark on in more detail in Section |H 



We will show that the even Pfaffian flow of an if 5 -tridiagonal matrix converges to 
a sequence of block diagonal matrices with the following shapes: 

(a) a 2 x 2 block for each real pair of eigenvalues (z, —z), 

(b) a 2 x 2 block for each imaginary pair of eigenvalues (z, —z), 

(c) a 4 x 4 block for each quadruple of complex eigenvalues ( Z ■ Z • Z • Z ) • 

In addition, for the tij flow, the blocks are sorted by the size of Re(2 2j '). In the appendix, 
we will give a simple proof that the SR-algorithm is given by the integer evaluation of 
the flow generated by the Pfaff lattice introduced in [1] (see also [16]). 
The present paper is organized as follows: 

In Section [2] we review the SR-factorization which is closely associated with the Pfaff 
lattice hierarchy of equations. We then show, in Section [3l that the odd members of 
the hierarchy of the Pfaff lattice (11.41) on symplectic matrices are trivial, while the even 
members of the hierarchy on if 5 -tridiagonal matrices are equivalent to the indefinite 
Toda hierarchy introduced in [I7J (Theorem 13. 2j) . This is similar to the case of the 
Toda lattice in the connection with the Kac-van Moerbeke system. In contrast with the 
symmetric Toda lattice hierarchy, the indefinite Toda hierarchy experiences blow ups, 
where some entries of the matrix approach infinity. Our result then implies that for 
generic if 5 -tridiagonal matrices the Pfaff lattice has blow ups. 

In Section HI we give a further discussion on the equivalence between the Pfaff 
lattice on if 5 -tridiagonal matrix and the indefinite Toda lattice defined in [17]. We 
begin by reviewing the indefinite Toda lattice, then show that the r-functions for the 
Pfaff lattice are the same as those for the indefinite Toda lattice (Theorem 14. II) . We then 
examine the families of skew-orthogonal polynomials appearing in the Pfaff lattice for an 
if 5 -tridiagonal matrix. We show that these polynomials are related to the orthogonal 
polynomials in the indefinite Toda lattice (Theorem 14.21) . Finally, we remark on the 
asymptotic behavior of the even Pfaff lattice flows of an if s -tridiagonal matrix. In 
particular we show that initial conditions with complex eigenvalues, or with complex 
eigenvectors, will result in a blow up. 



2. Background on SR-factorization 



The Pfaff lattice equation is intimately connected to the SR-factorization of an invertible 
matrix and so we collect here some relevant facts on this factorization. With the Lie 
subalgebra t given by (11.31) in the Lie algebra splitting sl(2n) = t © sp(n), we define G t 
to be the Lie group with Lie algebra t, 



( oixh 
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We also define the group G D Gf> to be the group of real invertible lower 2x2 block 
triangular matrices with free invertible diagonal blocks, that is, G is a parabolic subgroup 
of lower 2x2 block matrices of SL(2n, R). 

We define the Pfaffian of a skew-symmetric matrix m by the recursive formula 

2n 

pf(m) = ^(-l)^ +1 m uP f(my), 

where m--- is found by deleting the z-th and j-th rows and columns of m, where % is any 
row, and the initial value for the recursion is given by 

I m 12 \ 
pf = mu ■ 

V -mu o y 

Then, for example, with the skew symmetric matrix 

/ 777.12 m 13 m 14 \ 

m 23 m 2 4 
m 34 

V / 

we have 

pf(m) = m l2 m u - m n m 24: + m u m 22> . 

Throughout this paper, we leave blank the lower triangular part of skew-symmetric 
matrices. 

We will use the SR-factorization of g G SX(2n, R) introduced by Bunse-Gerstner 
(Theorem 3.8 in [6], see also [5]): 

Theorem 2.1 Let g G SL(2n, M) and M = gJg T . Then the factorization 
g = rs 

with r G G and s G Sp(n, R) exists if and only if the Pfaffian of the 2k x 2k upper left 
submatrix of M, denoted by M 2 k, does not vanish, i.e. 

pf (M 2fe ) ^ , 1 < k < n - 1 . 

As a consequence there is a dense set of matrices in SL(2n, R) for which this 
decomposition exists. If we restrict the factorization to using elements of G{ C G 
the Theorem becomes: 

Theorem 2.2 Let g G SL(2n, R) and M = gJg T '. T/ien i/ie factorization 
g = rs 

with r 6 G( and s G Sp{n, R) exists if and only if all the Pfaffians pf(M 2 fc) satisfy the 
positivity condition, 

pf (M 2fc ) > , 1 < jfe < n - 1 . 



This implies that there is an open set of matrices in SL(2n,M) for which this 
decomposition exists. 

Theorem 3.7 in |6j shows that the set of matrices in SL(2n, C) with SR-factorization 
is neither open nor dense in SL(2n, C). For this reason we will restrict our consideration 
to the real groups for now. 

In Remark 3.9 in [6] it is noted that for the choice of subgroups of SL(2n, M.) in 
that paper the decomposition g = rs is not unique, as one may pass a factor from 

G n Sp(n, R) = {diag 2 (A u A 2 , . . . , A n ) : det(A,-) = 1 Vj } , 

between r and s. In the case of SR factorization with R restricted to G% we find that 
factorizations are unique up to a factor from 

Gi n Sp{n, R) = {diag 2 (±/ 2 , ±/ 2 , . . . , ±/ 2 )} . 

One notes that G^ is in fact made up of 2 n connected components, and if we further 
restrict R to the component of G% containing the identity matrix we obtain unique 
factorizations. Equivalently we obtain unique factorizations by asking that the diagonal 
elements of r be positive. In contrast the group G is connected. 

The Pfaffians defined in Theorem 12.11 are also called the r-functions, which play 
the fundamental role in the Pfaff lattice hierarchy [HQS], i.e. for M = gJg T and M 2k , 
the 2k x 2k upper left submatrix of M, 

r 2k := pf (M 2k ) , fc = 0,l,...,n, (2.1) 

with t = 1. Derivatives of the r-functions for g(t) := exp {^^L^ tjL 3 Q j , with 
L(0) = Lq, generate the matrix entries of L(t) in the form (11. 5ft . and the solution 
of the Pfaff lattice hierarchy. For example, we have [HIS] (see Section [3]), 

a*(t) = a fe (0)^±SS35, hk(t) = JL i nT2fc(t) . (2 . 2) 
We also have the following Lemma: 

Lemma 2.1 Let g = rs be a 2n x 2n matrix where r £ G% and s G Sp(n, M). Then the 
2x2 block diagonal entries of r, diag 2 (r) = diag 2 (r 1 J 2 , . . . ,r n I 2 ), are expressed by the 
T-junctions, 



T2k 
T 2 k-2 



Proof. Since M = gJg T = rsJs T r T = rJr T , we have T 2k = pf(M 2 &) = nj=i r |' 
which leads to the formula. □ 

A necessary condition for the existence of the factorization is that T 2k = pf(M 2 fc) must 
be positive (Theorem 12.21) ; Lemma 12.11 and the convention that tq = 1 explain this 
condition. There are explicit algorithms for carrying out SR-factorization and they 
have many features in common with those developed for QR-factorization, we refer the 
reader to [9] and [6]. 



3. Pfaff lattice hierarchy 



In PP, Adler, Horozov, and van Moerbeke introduced the Pfaff lattice hierarchy to 
describe the partition functions of the Gaussian orthogonal and symplectic ensembles 
of random matrices (GOE and GSE), and to describe the evolution of the recursion 
relations of skew-orthogonal polynomials. The finite Pfaff lattice hierarchy is a set of 
Hamiltonian flows on sl{2n, R) with the Lie-Poisson structure induced by the splitting 
si(2n, R) = t © sp(n), and it is defined as follows (see also [IE]): For any X, Y e s[(2n) 
we first define the 7?.-bracket 

[X,Y]k = [KX,Y] + [X,KY], 

where 7Z is the R- matrix given by 7Z = — vr sp ). Then the Lie-Poisson bracket for 
any functions F and H on sV{2n) = sl(2n) is defined by 

{F, H} n (L) = (L, [VF, Vtf U , 

where (A, B) = tr(AB) and VF is defined by (X,VF) = f e F(L + eX)\ e=0 . The Pfaff 
lattice with respect to a Hamiltonian function H(L) is defined by 

^ = {H(L),L} n (L). 

In particular, if the Hamiltonian is Sp- invariant, one can write 

{H(L),L} n (L) = [Tr sp (VH(L)),L] . 

The traditional Hamiltonians giving the Pfaff lattice are ^-tr(L fc+1 ) [IE] . The Pfaff 
lattice hierarchy is then defined by 

— =[n 5p (L k ),L] k = l,...,2n-l. (3.1) 

Otk 

Thus the t^-flow of the Pfaff lattice hierarchy is associated to the Hamiltonian = 
^j-j-tr(L fc+1 ). Each flow is solved by the following factorization procedure: Factor 

exp (tV#(L(0))) = R(t)S(t) (3.2) 

with the initial matrix L(0) using SR-factorization with R in the connected component 
of Gi containing the identity and S G Spin), then the solution is given by 

L{t) = R^L^Rit) = S , (t)L(0)S'(t)- 1 . (3.3) 

With this in mind, one sees that as in the case of the QR-algorithm with the Toda lattice, 
the SR-algorithm is given by integer evaluations of the Pfaff flow with Hamiltonian 
H{L) = tr(Lln(L) — L) a result we will show in detail in the Appendix. 

In |16j . we showed that the Pfaff lattice hierarchy (13.11) is an integrable system in 
the Arnold-Liouville sense. Normalizing the matrix L of the form (11.51) by 

L = PLP- 1 

with the diagonal matrix P in the 2x2 block form, 
P = diag 2 J 2 , ail 2 , (aia 2 )h, ■■■ ■ 




L is a lower Hessenberg matrix with l's on the super diagonal. Then we showed that 
if L was not symplectic and all the eigenvalues are real and distinct, then as tj — ► oo, 
L(tj) converges to a 2 x 2 block upper triangular matrix such that the diagonal blocks 
are sorted by the size of zK 

Solutions of the Pfaff lattice hierarchy are generated by the r-functions (introduced 
as obstructions to the SR-factorization in the previous section). They are found by the 
following procedure: We first consider the factorization of g(t) := exp ^^j™! 1 tjL 3 
with L = L(0) (see (EOD ). 

g{t) = R(t)S(t) with R(t) G G h S(t) G Spin) . 
Then the skew-symmetric matrix M(t) = g(t)Jg T (t) becomes 

M(t) = R(t)JR T (t). 
Since R G Gf, we have 

diag 2 (i?) := diag (n h, r 2 h, • • • , r n I 2 ) ■ 
Then the r-functions r 2 k defined in (12.11) can be written by 

k 

r 2fc (t) = pf(M(t) 2fc ) = n^W 2 - 

3=1 

Then from (IQl . i.e. i2(t)L(t) = L(0)i?(t), we have 



a fc (t) = afc(O) -— = a fc (0)- v 



^fe(t) r 2 fc(t) 
which gives the a^'s in (12.21) (see also [U fT6]). 

Pfaff' flows of symplectic matrices 

Here we show that if L(0) is a symplectic matrix then it is a fixed point of the odd 
members of the Pfaff lattice hierarchy. To see this, one notes: 

Lemma 3.1 For L G sp(n), the odd power L 2 - 5-1 is also symplectic. 

Proof. Being symplectic is equivalent to JLJ = L T . Suppose that JL 2 - J_3 J = 
(L T ) 2,? 3 . Then we have 

JL 2 i-\J = JL 2j ~ 3 JJLJJLJ = (L T f- 3 L T L T = (L T ) 2j_1 

so that the lemma is true by induction. □ 

Therefore, 7r s a(L 2j_1 ) = L 2 ^ 1 for L G sp(n), hence the odd members of the Pfaff 
hierarchy become trivial, i.e. 

Note in particular that all bk in (12.21) vanish, which is consistent with the form L in 



Theorem Appendix A. 5, that is, the diagonal elements are all zero. 



{H j ,L}{L) = [ir m (VH j ),L). 



We note that this also happens for the Toda lattice hierarchy: Recall that the 
Toda lattice equation for a symmetric matrix is based on the Lie algebra splitting 
sl(n) = b ©so(n), where b is the set of upper triangular matrices. With the pairing 
(A,B) := tr(AB) for A, B e sl(n), we have st(n) = sl{n)* = b* © so{n)* where 
b* = so(n) ± = Sym(ra) and so(n)* = b 1 - = n. Here Sym(n) is the set of symmetric 
matrices, and n is the set of strictly upper triangular matrices. Then the Lie-Poisson 
bracket for the functions F, G on Sym(ra) is defined by 

{F, H} (L) = (L, [VF, VH]) for L e Sym(n) . 

The Toda lattice hierarchy for L e Sym(n) is then defined by 
dL 

With the Hamiltonian functions Hj(L) = ■^-tr(L J,+1 ) for j = 1 
equation for L becomes 
dL 

dT 3 

with 7r S0 (X) meaning the projection of the matrix X on the so(n) component. 

Now if we extend the Toda lattice equation for the general L e sl(n, M), an analogue 
of Lemma 13.1 1 is true with so(n) instead of Sp(n). Then, for L e so(n) (i.e. skew- 
symmetric), the odd members of the hierarchy become trivial. However the even powers 
of L are symmetric, and the even members of the generalized Toda hierarchy for L 
are equivalent to symmetric tridiagonal Toda lattices. This case includes the Kac-van 
Moerbeke system |14j : Consider for example the following In x 2n skew-symmetric 
tridiagonal matrix, 



n- 



(3.4) 

1, the differential 
(3.5) 
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Then the Kac-van Moerbeke hierarchy may be expressed by 



dL 



Of 



the first member gives 



2.7 



[7r S0 (L 2j ), L], where 



dt 2 



a 



k+l)i 



k = 1, 



2n- 1 



with «o = «2n = 0. We then note that the square L 2 is a symmetric matrix, 



L 




(3.6) 



where Tw, for i = 1,2, are n x n symmetric tridiagonal matrices given by 



I bf 



r/ 
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n-1 






*n-l 
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with a 



(i> 



(2) 



and b 



(2) 



a 2fc-l a 2fc- 



«fc — a2fc-ltt2fc, Ofc — — a 2k-2 ~ a 2k-H a k — a 2k°i2k+li "k 

Then one can show that each gives the symmetric Toda lattice, that is, the equation 
dL 

[ir S0 (L i), L] splits into two Toda lattices, 



Of 



Of 



2j 



[ir„(T&y,T®] for £ = 1, 2 . 



The equations for T"W are connected by the Miura-type transformation, with the 
functions (a^\b^) } through the Kac-van Moerbeke variables (see [T3]). 

We now show that the even Pfaff flows on a if s -tridaigonal matrix have a similar 



structure. 



3.2. Even Pfaff flows of symplectic matrices 

We first show that the if 5 -tridiagonal form ( 11. 6ft is invariant under the even members of 
the Pfaff hierarchy. Then we show that the Pfaff flow for the matrix L in the form (11.61) 
is related to the indefinite Toda lattice defined in jTT] . The indefinite Toda lattice is a 
continuous version of the HR algorithm (TO] , and is defined as follows (see also Section 
H]for some details): Let C be the n x n diagonal matrix given by, 

C = diag(ci,c 2 , . . . ,c n ) , 

where c^'s are the elements appearing in the diagonal blocks of (11.61) . and they are 



invariant under the Pfaff flow 
triadiagonal matrix ( 1A.2I) . i.e. 

( C\d\ c\d\ 

c 2 a l C2OI2 

T " 



sec 



We also define T by T := CT with the 
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Cn— idn— 1 C n — l<2n— 1 
Cn^n— 1 C n d n 

Then the indefinite Toda lattice hierarchy is defined by 
df 



(3.7) 



where Ba 



dt 2j 
[f j }+ 



(3.8) 



[r j ]_. Here the [T J ']± represent the projections on the upper (+) 



and lower (— ) triangular parts of the matrix T- 7 (sometimes we write T± = [T^}±] 



With L in the if s -tridiagonal form (11.61) . we first note 

(3.9) 




from which we have, for any even power, 

L 2j = T* ® ( I ° ) + (fif ® 






Remark 3.1 VFe remark here that equation $3. 9\) is similar to $3. 6\) . this structure is 
fundamental and implies that our problem will have a similar structure to that of the 
Kac-van Moerbeke lattice. 

Then using the projection (11.21) . we obtain 

(3.10) 

with Bj — — T 3 _. With (13.101) . one can easily show that the if ^-tridiagonal form is 
invariant under the even members of the Pfaff lattice hierarchy. It is also immediate to 
see from ( 13. lOf) that the 2m-th member of the Pfaff lattice hierarchy (11.41) in terms of L 2 
is equivalent to the indefinite Toda lattice (13.81) . Thus we have the following theorem: 

Theorem 3.2 The even Pfaff flow for L e sp(n, M.) in the H -tridiagonal matrix form 

is equivalent to the indefinite Toda flow for T = CT with the symmetric tridiagonal 
matrix T of liA.ty) and C := diag(ci, c 2 , . . . , c n ), 



^ = [B m ,f], (3.11) 

where B m = T™ -f™. 

If the Cfc = 1 for all k, then T = T is symmetric and B m = 7r S0 (T m ). Theorem 
13.21 then says that the 2m-th Pfaff flow is just the traditional symmetric tridiagonal 
m-th flow of the Toda lattice hierarchy. In terms of the SR-algorithm in Proposition 



Appendix A. 4 , we obtain as a Corollary of Theorem 13.21 that the even-integer iterates, 
L 2 k of the SR-algorithm are equivalent to the Tj, iterates of the HR-algorithm. This 
result is shown directly in |4j, the key is the structure of equation (I3.9p . With the 

permuted J matrix used in [4J, i.e. J = ^ ^ J | ®/ n , the right hand side of relation 

(13. 9p takes the form 

f 




4. Pfaff lattice vs indefinite Toda lattice 



Here we give a further discussion on the equivalence between the Pfaff lattice hierarchy 
on i? 5 -tridiagonal matrix and the indefinite Toda lattice hierarchy introduced in (IT) . 

4-1. The indefinite Toda lattice 

Let us begin with a brief description of the indefinite Toda lattice defined by (13.81) for 
a matrix T = CT with a symmetric matrix T of the form ( 1A.2I) . i.e. 

^ = [B,f] with B = f + -f.. 

The solution T(t) with the initial matrix T can be solved by the HR-factorization, 

g = exp(tTo) = r(t) h(t) , 

where r is a lower triangular matrix and h satisfies hCh T = C (recall that if C = I n , 
h G SO(n)). It is then easy to show that the solution T(t) is obtained by 

T(t) =r(t)- 1 f r(t) = h{t)f Q h(t)- x . (4.1) 

We note that the entries a k in T are expressed in terms of the diagonal elements of r, 
say diag(r) = diag(ri, . . . , r n ), 

a k \t) = afc(O) 7— . 

r k {t) 

To find those diagonal entries from g, one considers the following matrix called the 
moment matrix, 

M Toda := ~gCg T = rCr T . (4.2) 

Then the matrix r can be found by the Cholesky factorization method. Now the r- 
functions are defined by 

Tfe Toda = det (M fc Toda ) , 

where M^ oda is the k x k upper left submatrix of M Toda . With (14.21) . we have 
r Toda _ n . =l C j r| , which gives c k r\ = r k /r k _i. Then the entries a k of T can be 
expressed in terms of the r-functions, 

As will be shown in the next section (see Theorem 14.11 below), those are the same as 
the formulae for the a k s in (12.21) given in terms of the r-functions of the Pfaff lattice. 

Let us also discuss the orthogonal functions appearing in the indefinite Toda lattice. 
First we note that the Lax form (13.81) is given by the compatibility of the equations, 

T$ = $L> and —— = M , 

dt 2j 

where D = diag(Ai, . . . , A n ) with the eigenvalues X k , and $ = (4>i(\j))i<ij< n the 
eigenmatrix. As the orthogonality condition of the eigenvectors with a normalization, 
we have 



Lemma 4.1 = C . 

Proof. We consider the following relation which is equivalent to the orthogonality 
relation, 

§ T c- l <$> = c- 1 . 

(This is sometimes called the second orthogonality relation.) To show this, we note 

A i 0(A i ) T C-V(A i ) = ^fC^f^) = 0(A,) T T0(A,) 

= (T0(A l )) T C- 1 0(A,) = A J 0(A J ) T C- 1 0(A,) . 

Since we assume \ 7^ Xj, the matrix $ T C~ 1 $ is diagonal. We also note that this matrix 
is invariant under the flow, i.e. d(& r C~ 1 &)/dt = 0. Normalizing the diagonal matrix 
gives the result. □ 

We then define the inner product for functions /(A) and g(X), 

n 

(/^)Toda := f( X k)g( X k)c k , (4.4) 
fc=l 

which defines a discrete measure dfi(X) = Ylk=i c k$(. x ~ ^fe) dz. Lemma I4TT1 implies 

(0j0j)Toda = CiSij . 

With the orthogonality relation, the moment matrix M Toda can be expressed by 

M Toda = ^B C< = ( (0 O O e2t2)Toda) ^ ( (45) 

where $0 is the initial eigenmatrix with the eigenvector 0° = 0° (A) = <pi(X, 0). 

In terms of the initial eigenmatrix $0, we can find an explicit form of $(t) for the 
indefinite Toda lattice hierarchy with t = (£2, t&, ■ ■ ■ , t2 n ) (see Theorem 1 in [17] and 
Theorem 2 in |!5j): 

Lemma 4.2 The eigenvector 0(A,t) = (0i(A,t), . . . ,0 n (A,t)) T can be expressed as 

mii(t) •■■ mi, fc _i(t) 0?(A) 



0fc(A,t) 



c V 2 e 5(A,t) 



m 2 ,i(t) •■■ m 2 , fc -i(t) 0°(A) 



(4.6) 



[r Toda (t)r Toda (t) ]^ 

m M (t) •■■ m fcjfc _i(t) 0fe_x(A) 

where rriij(t) := (0°</>°e 2 ^ A,t ))Toda £(A, t) = Ylk=i ^hk, and the T-junctions are 
given by 

r fc Toda = det (M fc Toda ) = \{m^) x <^< k \ . (4.7) 

Proof. With the factorization e s(To,t) = r(t)h{t), we have T(t) = /i(t)f /i(t)- 1 (see 
(14. ip ). Then from Lemma [4.11 we obtain 

$(t) = /i(t)$ . 

Now from the factorization with T$ = QD, we have 

$(t) = r(t)" 1 e« fo>t) $ = r(t)' 1 $ e* (D,t) , 



which implies 

4>i(\ t) G Span (0°(A)e^ A ' t ), . . . , $(\)e^ } , i = 1, . . . , n . (4.8) 

Then using the Gram-Schmidt orthogonalization method with the inner product (I4.4p . 
we obtain the result. □ 

If we set 0?(A) = 1 and consider the semi-infinite lattice with n = oo, fc (A,t) can 
be expressed in the following elegant form in terms of only the r-functions: 

Proposition 4.1 The orthogonal eig en functions (j) k {\) can be expressed in terms of t- 
functions, 



0fc(A,t) 



_Toda 



[r, Toda (t) « a (t)] 



1/2 'k-1 



t-itA" 1 ])^- 1 



wfcere £(A, t) = A fc * 2fc and r/ oda (t - ifA- 1 ]) 



Toda 
j 



t-2n 



1 



2nA" ' 



•)• 



Proof. This proposition appears in [2] for the case c k = 1. In this more general 
case the proof follows from formula (14.61) with £(A,t) = YlT=i A fc ^2fc: First we note 
that using (14.81) . one can replace 0°(A) in (14 .6p by A J_1 with 0?(A) = 1, which 



gives mij(t) 



(A 



i+i-2„2f(A,t) 



)Toda- This implies that mjj(t) = d %+ ^ 2 m 11 (t) / dt\ 



hence TjF oda (t) is given by the Hankel determinant form, r^J 



Toda 



\(h i+j-l)l<i,j<k\ 



with h 



i+j-i 



rriij. Since hj(t) is a linear combination of the exponential function 



Ei(t) = e 2 ^^ an d E^t - ifA" 1 ]) = (1 - ¥)Ei(t), we have 



Then it is easy to see that 



_Toda 



1 



[A- 1 ] A 



-h j+1 (t) 



h 2 (t) 



h 2 (t) 
h 3 (t) 



h k +i(t) 
h k+2 {t) 

-Toda 



h k (t) 

h 2 k-i(t) A 



1 

A 



fe-i 



h 2k (t) A A 



h k {t) 
h k+ i(t) 

which gives the determinant in (14. 7p in terms of r^_ 1 . 

For a finite n, one can see that r^ oda (t — |[A _1 ])A n /rJ oda (t) is proportional to the 
polynomial rjfe=i(A — A^), and we have n+ i(A) = 0, which is just the characteristic 
polynomial of L. □ 

With Proposition I4.1[ we will show a further example of the close relation between 
the orthogonal functions in the indefinite Toda lattice and the skew-orthogonal functions 
in the Pfaff lattice for the symplectic matrix of (11.61) . 



4-2. The t -functions of the Pfaff and Toda lattices 

We show here that the r-functions which generate the solutions of the Pfaff lattice 
equations are equivalent to the r-functions which generate the solutions of the indefinite 



Toda lattice hierarchy. Let us recall that the r-functions of the Pfaff lattice are defined 
by fl2U), i-e. 

r 2k = pf (M 2fc ) , 

where M 2k is the 2k x 2k upper left submatrix of 2n x 2n skew-symmetric matrix M 
given by 

M:=gJg T with g = e c(L o' t} , 

where £(Lo,t) = Y^k=i ^2kLl k with the initial matrix L(0) = Lq. Then our goal is to 
show the following Theorem: 

Theorem 4.1 The r functions of the even-flows of the Pfaff lattice hierarchy with the 
initial matrix L(0) in the H -tridiagonal form of U.6\) are related to the r-functions of 
the indefinite Toda lattice hierarchy by 

1 Toda 
T 2k = T k 

Proof. We first note from (13.91) that the entries of g = e^ L °'^ are expressed by those 
of g = e€( f o,*), 

„ ( 1 \ „ji / 

y y \ooy y \oi 

Since T = CT with a symmetric matrix T , C -1 ^ is symmetric, i.e. C~ l g = g T C~ x . 
Then for H = C^M™* = C~ l gCg T = g T g T and M = gJg T , we have 

Let i^fc be the k x k upper left submatrix of H, and M 2k be the submatrix defined 
above. Then with the permutation matrix 

P = [ei, e 3 , . . . , e 2 fc_i, e 2 , e 4 , • • • , e 2fc ] , 

where e^'s are the standard basis (column) vectors on M. 2k , one can transform M 2k to 

"«*-(-%: I) pTjp {tl)- ^ 

We recall the Pfaffian identity for a skew-symmetric matrix B, 

pf(A T BA) = det(A)pf(B) , 
which with (I4.10p implies that 

det(P)pf(M 2fc ) = det(P)det ( ^ °* j pf(J) . 

We finish the proof by noting that P is invertible and pf( J) = 1 giving 

l>f(-U 2 ,)=<H.( H " M =det(PT fc ) = ^^det(M fc Toda ). 
fc 4 / ci ■ • • c fc 



4-3. Skew- orthogonal and orthogonal functions 

We show here that the eigenvectors of the matrix L in the form ( 11.61) define a family 
of skew-orthogonal functions, and then discuss an explicit relation between those skew- 
orthogonal functions and the orthogonal functions appearing in the indefinite Toda 
lattice, i.e. those in Proposition 14. 1L 

As in the case of the Toda lattice, the Pfaff lattice hierarchy (13. ip is given by the 
compatibility of the linear equations, 

L^ = ^A and J^- = vr sp (L 2j ) , (4.11) 

dt 2j 

where L is symplectic and in the form of (11.61) . and A is the eigenvalue matrix, 

with D = diag(^i, . . . , z n ). (Recall that the eigenvalues of L consist of the pairs (zk, —z k ) 
for k — 1, . . . , n, and we assume here that they are all distinct.) Then we have: 

Lemma 4.3 The eigenmatrix \I/ satisfies the skew- orthogonal relation, 

# T = icj, 

with a diagonal matrix with nonzero constants K k 's in the form, 
JC = diag 2 («i/ 2 , ^h, • • • , Knh) ■ 
Proof. From (14.111) with B := 7r sp (L 2j ) and t = t^j, we have 
d 

(m T jm\ = m T B T jm + m T jBm 

at 

= m T (B T J + JB)m = 0, 

where we have used B e sp(n). 

Now we note the following equation for the eigenvectors ip(z k ), 

Zjtpizk) 7 'Jip(zj) = il){z k ) T JLip{z j ) 

= -ij(z k ) T L T J^( Zj ) (y JL + L T J = 0) 

= -{mz k )Y 

= -z k ^(z k ) 7 'Jip(zj) . 

This implies 

(z k + Zj)tjj(z k ) 7 'Jtp(zj) = . 

Noting the order of the eigenvalues in A = diag(— zi, zi, . . . , — z n ,z n ), so that \1> T J\I/, 
which is skew-symmetric, is a multiple of J of the form /C with nonzero constant K k s. 
□ 

From Lemma 14. 3} we have 

tf£ _1 Jtf T = J, (4.12) 



which implies the skew-orthogonality relation for the functions ipk{ z ) in the eigenvector 
ip(z) = (ipo(z), ipi(z), . . .,-02n-i(^)) of L, 

(^•,« = (^ + i,^ +1 ) = o, for 

Here the inner product (/, g) is defined by 

n 

(/, g) ■= Yl lf(- z Mz k ) - f{z k )g{-z k )\ k- 1 . (4.14) 

k=l 

Now we claim: 

Theorem 4.2 The skew-orthogonal eigenf unctions ipk{z,t) satisfying l[4-13{ ) can be 

given by 

1p2h{z,t) = 0fc+i(A,t) 

with A = z 2 , 

ifak+i(z,t) = c^z <j) h+x {\t) 

where 0fc(A, t) 's are the orthogonal eigenf unctions appearing in the indefinite Toda lattice 
(see Proposition 4-1), and the measure of the inner product (•, •) is given by n k = 2z k c k ~ l ■ 

Before proving the Theorem, we recall the following Proposition which gives the 
skew-orthogonal functions appearing in the semi-infinite Pfaff lattice: 

Proposition 4.2 The eigenvector ip(z,t) = (ipo(z, t), ip\{z, t), . . .) T for the matrix L 
of U.5\) with n = oo can be expressed in terms of r-functions, 

[r 2fc (t)r 2fc+2 (tO] 1/2 



^2k(z,t) = — ——j-r^t ~ [Z X ])z 2fc 



where £(z,t) = Y^kLi^k and ^(t - [z^ 1 ]) = T 2 k(h - ~, t 2 - ^, . . .). 

This Proposition appears as Theorem 3.2 in [3] for the case = 1. In this more general 
case the additional factor of a c^ +1 is present to ensure that the recursion relation L has 
L 2 k+i,2k — Cfc+i- This entry of the Pfaff variable is a Casimir of the Pfaff lattice equations 
and so is fixed. In the definition of the skew-orthogonal functions it corresponds to a 
choice of the ratio between leading coefficients of the polynomial parts of ip2k and ip2k+i- 
Now we prove Theorem 14.21 

Proof. First we show that ^>&(z, t)'s in the Theorem are the same as those in 
Proposition 14.21 when the matrix L is in an if 5 -tridiagonal form (jl.6p . To show this, 
we recall Theorem 14.11 i.e. 

Cl • • • C k 



Since the odd flows are trivial, we have t = (t 2 , t^, . . . , t^n)- Then substituting this 
relation into ip k (z,t) in Proposition 14.21 it is straightforward to show the equations in 
the Theorem. 

Now we show the skew-orthogonality relation (14.31) : We only check the case 
{4>2j, i>2k+i) — djk, and the others are trivial with the form of the inner product. 

n 

(i>2j,i>2k+l) = ^2 [^•(-^V'aJfc+l^) - 1p2j{Zi)lp2k+l{-Zi)] 7^ 
i=l lZi 



= C k+1 ^2 0i+l(Ai)0fe+l(Ai)Cj = S jk , 
i=l 

where we have used the orthogonal relation $C$ T = C in Lemma [4. II □ 

We also note that the moment matrix M(t) for the Pfaff lattice has the similar 
form as the M Toda (t) for the indefinite Toda lattice given in (14.51) . 

M(t) = e«*> Je m*) T = y e^ A2 V *- 1 ;$- T e«( A, ' t ) 

\ J / 0<i,j<2n-l 

Comparing with (14. 9p . the entries m^- := (ip®, ipVe^* 2 '^) are given by m^- = — rriji with 
m 2j -i,2k = -m 2k ,2j-i for j < k , 

and zero for all other cases, i.e. m 2 j-i t 2k~i — , m>2j,2k = 0. With this structure of the 
moment matrix, one can give a direct proof of Theorem 14.21 without Proposition 14. 2\ 
however this might be a less elegant approach. 

Remark 4.3 The skew-inner product of is closely related to the skew-inner 

product for GSE- Pfaff lattice, which is defined as 

n 

{fi 9) GSE := ^2{f,9Kzzk)c k , 

k=l 

where {f,g}(z) = f'(z)g(z) — g'(z)f(z) with f'(z) = df(z)/dz. This inner product may 
be obtained by the following inner product in the limit z 2k -\ —> Z2k (see JJE/), i.e. 

,, v r f(Z2k-l)g(z2k) ~ f{z2k)g{z2k-l) 

{f, 9)gse = am > c fc , 

fc=1 Z 2k — z 2k-l 

In the case of ^.14\ ), we instead take the limit z 2 k-i —* —Z2k- 

Also in the case of the continuous measure, the inner product ( [^. 1$ can be written 

by 

(f,g)= [ [f{-z)g{z)-f{z)g{-z)\ c(z 2 )dz, 

where S = M. + U zR + oriented from ioo to and then to 00. To reduce this expression 
back to the discrete case we take 



c(z 2 ) dz\^ = ^2 $(z 2 ~ z l) dz = ^2 $(( z ~ z k)(z + Zk))dz^ 
k=i E k=i s 

n r 1 1 1 n 1 

= —5{z-z k ) + —5{z + z k ) dz = ^ — 5(z - z k ) dz 
z,z k £z k x 



fc=l 



k=l 2zk 



with Zk G S, and where we used the property 5(az) = m<£(z) /or a G R and 
5(az) = ihjfi(z) f or a G zR. The last equality is because we have restricted to the 
positive and positive imaginary square roots with Zk G E. The result is that this discrete 
inner product will agree with fl^.i^p . 

Asymptotics of the even Pfaff flows 

In this final section, we mention the asymptotic behavior of the even Pfaff lattice based 
on the results of the indefinite Toda lattice discussed in [T71 [18] : Theorem 13.21 implies 
that the Pfaff flow on iT 5 -tridiagonal matrix Lq with Ck = ±1 is equivalent to the 
indefinite Toda flow on T = CT given by (13. 7p . This system was studied in detail 
in [T7] and [12]. It is a version of the Toda lattice which uses HR-factorization in 
place of QR-factorization, see [10J. The goal of HR-factorization is to write an element 
g G SL(n, R) as g = rh where r is lower triangular and h satisfies 

hCh T = C . 

For g(t2j) := exp(t2jT(0y), the HR-factorization g(t 2 j) = r(t2j)h(t2j) gives the solution 
of the t2j-flow of the Pfaff lattice, 

f(t 2j ) = r-^f (O)r(fjy) = h(t 2j )f (O)^- 1 ^) . 

Thus the indefinite Toda lattice is a continuous version of the HR-algorithm [TO] . 

The indefinite Toda flow may experience a blow up, where some of the entries reach 
infinity in finite time [18]. A blow up occurs when one of the r functions becomes 0. In 
our case this is precisely when one of the T2k = 0. The initial conditions for a blow-up 
are characterized by 

Theorem 4.4 (Theorem 3 in [17J) IfT(0) possess non-real eigenvalues or non-real 
eigenvectors while r 2 fc(0) ^ 0, then T{t) blows up to infinity in finite time. 

The second possibility occurs if the Ck do not all have the same sign. The case when the 
eigenvalues of T are real and the Ck do not all have the same sign is studied in Section 4 
of [T5] . It is shown that every flow (considered in both directions t 2 j — > ±oo) contains a 
blow up. We may then compactify the flows by adding infinite points representing the 
blow ups. It is then shown that the fixed points of the compactified flows are diagonal 
matrices. In other words — > and dk — > CkZ^, k s as t 2 j — > oo for some permutation 
a G S n , the symmetric group of order n. When the eigenvalues of T are real we may 
count the total number of blow ups on the flow. If the number of changes of sign in the 
sequence Ck is m then the number of blow ups is m(n — m) [18J. 

There is a condition under which the indefinite Toda flow, with Ck = — 1 for some k, 
does not contain a blow up in the positive t direction. A lower triangular matrix is called 
lower triangular totally positive if all its non-trivial minors are positive. It was shown 
in [12] that if the lower component of the LU-factorization of the eigenvector matrix of 
T is lower triangular totally positive then there are no blow ups on the indefinite Toda 
flow for t > 0. 



If L has only real and imaginary eigenvalues, then T has real eigenvalues. By 
Theorem 13.21 we find that the flow may be continued through a blow up, if it occurs, 
and that this compactified flow converges: 



where the eigenvalues of L are {±z m : m = 1, . . . , n}. 

If the Cfc = 1 for all k, then T is symmetric and no blow up is possible. We will 
restrict to this case from now on. We appeal to the abundant literature on the Toda 
lattice for the results needed [El [T71 fl8] . In particular, under mild assumptions, as 
t 2 j — > ±oo, T(t 2 j) given by (13.111) converges to a diagonal matrix. As a result we see 
that CLkfoj) — > 0. 

In terms of the Pfaff lattice this is showing that £(£2) converges to a 2 x 2 block 
diagonal matrix as t 2 —>■ 00 where the j-th diagonal block has the form 



The blocks will be sorted by the size of z 2 so that the blocks with pairs of imaginary 
eigenvalues will appear in the lower right corner while those with pairs of real eigenvalues 
will appear in the upper left corner. For the inflow, the eigenvalues will again be such a 
pair but will now be sorted by the size of z 4 , mixing the blocks with real and imaginary 
pairs. 

While the end result is that the Pfaff lattice restricted to symplectic lower 
Hessenberg matrices is equivalent in every sense to the indefinite Toda lattice flows, 
this is not a result which is immediately obvious from the Pfaff lattice itself. One sees 
that even in the case when the Ck = 1 for all k, the long time dynamics of the Pfaff 
lattice flows of symplectic lower Hessenberg matrices is significantly different from those 
of non symplectic lower Hessenberg matrices; for example the the symplectic matrices 
are in fact fixed points of the flow. This difference alone merited a separate work on 
these cases and produced some interesting and new connections to generalizations of 
the Toda lattice. 

Appendix: SR-algorithm 

The Toda lattice may be viewed as a continuous version of the QR-algorithm for 
diagonalizing a symmetric tridiagonal matrix. In the same way the Pfaff lattice is a 
continuous version of the SR-algorithm on an if 5 -tridiagonal matrix. Here we collect 
some pertinent details about the SR-algorithm. These results have been covered from 
a variety of points of view, a partial list of references is [U [71 EE [191 1201 [ZD [221 [23] . 

Here we consider the SR-algorithm for a symplectic matrix which is symplectically 
similar to a lower Hessenberg matrix. Recall that any matrix can be reduced to a similar 
Hessenberg matrix by Householder's method, a fact used in giving an efficient version of 





the QR- algorithm. However, note that the Householder's method does not, in general, 
give a symplectic conjugation. 

We will need Theorem 3.4 from [7j: 

Theorem Appendix A. 5 Let h 6 Sp(n, M) and Si be a row vector in R 2n . Then there 
exists a symplectic transformation S such that ShS~ x is a lower Hessenberg symplectic 
matrix iff K{h, si) has an SR-factorization K(h, si) = RS where R is a lower triangular 
matrix and Si is the first row of S G Sp(n, M). In addition, if this factorization exists, 
then 



Lq := ShS = R ChR 
is in the lower Hessenberg form, 
( 



ci 
d 1 



ai 




o 2 




ai 


c 2 
d 2 




o 2 










o 2 


o 2 




c n 

d n 



V 

that is, Lq has a 2 x 2 block tridiagonal form Lq = (/ 



(A.1) 



/ 



i,j)x<i,j<n with 2x2 block matrices 



lij having lij = 2 for \i 




k,k+l 




where Ck = ±1. 

A dense set of symplectic matrices h may be placed in this form. We call a matrix in 
the form (jl.6p an "if 5 -tridiagonal" matrix, which plays the similar role as a tridiagonal 
matrix in the case of the symmetric matrices. 

Furthermore, at the expense of excluding a large set of symplectic matrices, we may 



refine Theorem Appendix A. 5 



Theorem Appendix A. 6 Let h 6 Sp(n, M), and si be a row vector in M. 2n . Then there 
exists a symplectic transformation S such that ShS -1 is a lower Hessenberg matrix iff 
K(h, si) has an SR-factorization K(h, si) = RS with R e G%, and with s\ equal to the 
first row of S . In addition, if this factorization exists, then 

Lq := ShS 1 = R ^ChR 

is in an H s -tridiagonal form U.6]) with = 1 for all k. 

We now show that if h has quadruples of complex eigenvalues then it will not be 
similar to a matrix in the if ^-tridiagonal form of (11.6j) with Ck = 1 for all k: 

Proposition Appendix A. 3 Let L be a2nx2n matrix in the form of U.6\) with Ck = 1 
for all k. The eigenvalues of L come in real or imaginary pairs without multiplicity. 



Proof. One checks that for L in this form, we have 
L 2 = T®I 2 

where T is a symmetric tridiagonal n x n matrix given by 
/ dt ai • ■ ■ \ 

CL\ G?2 0*2 ' ' ' 

• • • d n _! a n _i 
V • ■ • a n _i d„ / 



(A.2) 



From this one computes that 

F(X) = det (L 2 - XI 2n ) = [det (T - XI n )} 2 . 

As T is symmetric, F(\) only has n real roots each with multiplicity 2, therefore the 
eigenvalues of L are square roots of real numbers, and so are only real or imaginary. □ 



There is an algorithm for carrying out the transformations of Theorems Appendix A. 5 



and Appendix A.6 , see [9]. 

This procedure is the analogous operation to the QR-algorithm step of transforming 
a symmetric matrix to a tridiagonal form by the Householder's method. There is a large 
time savings in carrying out this transformation first as an SR-factorization of a matrix 
in the form of L Q only requires 0{n) operations. 

The SR-algorithm is defined as an iteration with initial matrix Lq a symplectic 
matrix and recursion given by factoring L^-i = RkSk using the SR-factorization then 
taking Lk ■= SkRk = SkLk-iS^ 1 = R^Lk-iRk- As each step is a similarity transform 
of the previous step by both a matrix in the symplectic group and a matrix in Gt, Lk 
is still a lower Hessenberg matrix which is also symplectic, and therefore is still in the 
if s -tridiagonal form of L above. 

If the algorithm is successful the Lk approaches a family of block diagonal matrices 
with blocks of the form: 

(a) 2x2 blocks containing two real eigenvalues (z, —z), 

(b) 2x2 blocks containing two imaginary eigenvalues, (z, —z), 

(c) 4x4 blocks containing a quadruple of complex eigenvalues ( 

In addition the blocks are sorted by the size of hi \z\. If there are complex eigenvalues the 
sequence does not converge to a fixed matrix, rather it approaches the sorted diagonal 
shape. If the Cj = 1 for all j, then Lk converges to a block diagonal matrix with just 
2x2 blocks. The rate of convergence is 0(k 3 ) for a dense set of L . In practice one 
runs the algorithm until sup^- |aj(A;)| at the kth step is less than some fixed e tolerance. 
The algorithm also works on the initial matrix h (i.e. without the change to a lower 
Hessenberg L ). 

There is a substantial literature on improvements to this basic algorithm using 
implicit SR-factorization steps on certain matrix functions of Lk rather than just Lk 
(see e.g. [9] and references therein). 



We now show that the SR-algorithm is directly equivalent to the Pfaff lattice flows 
with a non-traditional Hamiltonian. In light of the connection to the indefinite Toda 
lattice and HR-algorithms this is not a surprising fact and the proof is appropriately 
close to that for the analogous fact in the Toda cases. 

Proposition Appendix A. 4 Let Lq e sp(n) in the form of ( f 1.5{) . Then the SR- 
algorithm is equal to the integer evaluations of the Pfaff lattice flow with respect to 
Hamiltonian H(L) = tr (Lln(L) — L) with L(0) = Lq. 

Proof. Recall that the SR-algorithm for the initial matrix Lq is given by 

Lk-i = RkSk and L k := SkRk ■ 

One can see that 

LC C C C T C — 1 C — 1 C — 1 G~ 1 

k — Ok^k-l " " " O20i-Lo0 1 D 2 ■ ■ ■ D k _ 1 D k . 

While the Pfaff flows arise from (13.31) . 

L(t) = S(t)L S(t)- 1 

where 

exp (t ln(L )) =L\ = R(t)S(t) . 

We want to show that S(k) = S^S^-i • • ■ S 2 Si. We prove this by induction: First we 
check that 5(1) = Si, from which L\ = L(l). Next we make the inductive hypothesis 
that S(k — 1) = Sfc-iSfc_2 ■ ■ ■ S 2 Si and then consider 

R(k)S{k) = exp(Hn(L )) = L k = L^Lq 
= R(k-l)S(k-l)L 
= R(k - l)L k _ 1 S(k - 1) 
= R(k-l)R k S k S(k-l) . 

By uniqueness of SR-factorizations for R in the identity component of G%, we see that 
S(k) — S k S(k — 1) = SfcSfc-i • • • S 2 Si. □ 
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